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Abstract. The ISM, powered by SNe, is turbulent and permeated by a magnetic 
field (with a mean and a turbulent component). It constitutes a frothy medium that is 
mostly out of equilibrium and is ram pressure dominated on most of the temperature 
ranges, except for T< 200 K and T > lO^K, where magnetic and thermal pressures 
dominate, respectively. Such lack of equilibrium is also imposed by the feedback of the 
radiative processes into the ISM flow. Many models of the ISM or isolated phenomena, 
such as bubbles, superbubbles, clouds evolution, etc., take for granted that the flow 
is in the so-called coUisional ionization equilibrium (CIE). However, recombination 
time scales of most of the ions below 10^ K are longer than the cooling time scale. 
This implies that the recombination lags behind and the plasma is overionized while 
it cools. As a consequence cooling deviates from CIE. This has severe implications 
on the evolution of the ISM flow and its ionization structure. Here, besides reviewing 
several models of the ISM, including bubbles and superbubbles, the validity of the CIE 
approximation is discussed, and a presentation of recent developments in modeling the 
ISM by taking into account the time-dependent ionization structure of the flow in a 
full-blown numerical 3D high resolution simulation is presented. 



1. Introduction 

In the last decade due to the substantial increase of computing power, sophisticated 
models of the interstellar medium (ISM), including the disk-halo interaction, bubble 
and superbubble evolution, molecular clouds evolution and fragmentation, formation 
of shock compressed layers, etc., in a turbulent supernova-driven medium, have been 
developed. In general these models included cooling and heating, which determine the 
physical state of the gas. Interstellar cooling can be the result of line and continuum 
plasma emission processes, as well as of adiabatic expansion of over-pressured gas. 
The importance of the former depends on the amount of atoms and ions present in the 
flow, whereas the latter is related to thermodynamical processes. For optically thin in- 
terstellar plasmas, frequently the assumption of coUisional ionization equilibrium (CIE) 
is used, according to which the coUisional excitation of the gas is followed by photon 
emission, with the number of ionizations being equal to the number of (dielectronic or 
radiative) recombinations. However, CIE violates detailed balancing, since colUsional 
excitation involves three particles, an atom or ion, an electron to collide with, and a 
second electron which is ejected, whereas in radiative recombination, the third particle 
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is a photon, leaving the system. Thus, strictly speaking, CIE can never be maintained 
in an evolving plasma, although it might be a fair approximation, especially in hot 
(T> 10^ K) environments. For lower temperatures recombinations are not synchro- 
nized with the cooling and therefore, deviations from CIE inevitably occur (see e.g., 
Kafatos|1973t |Shapiro, & Moore|1976l|Schmutzler & Tschamuter|1993| ). 



The CIE assumption in all ranges of temperatures, or at least in part of it, in the 
models discussed below implies the use of a specific cooling function, for given abun- 
dances, in the computational domain and at each time step. Obviously, the history of 
the ionization structure, which depends on the thermodynamic path of the plasma, will 
be lost, in contrast to a full non-equilibrium ionization (NEI) structure giving rise to a 
time-dependent cooling function, which is quite different than that obtained under CIE 
conditions. 

The structure of this paper is the following: In Section 2 a review on ISM mod- 
elling including superbubbles and disk-halo interaction, is presented. Sections 3 and 4 
deal with the collisional ionization conditions adopted in the models, and compare CIE 
with NEI results, respectively. As cooling is a process that depends on the history of the 
plasma, we discuss in Section 5 the cooling functions and emission spectra of gas in the 
latest simulations of the ISM. The paper is closed with a final remarks and conclusions 
section. 



2. From the 3-Phase Model to Present day ISM and Superbubbles Models 



Theoretical studies during the last three decades have culminated in the supernova reg- 



ulated Interstellar medium models of Cox & Smith ( 1974 1 and in the widely accepted 
"standard-model" ( McKee & Ostriker|1977 ), in which the gas is distributed into three 
phases in global pressure equilibrium, a cold and warm neutral phase (CNM and WNM, 
respectively), a warm ionized (WIM) and a hot intercloud (HIM) medium. There is 
global mass balance by evaporation, ionization and condensation, and energy balance 
between supernova (SN) energy injection and radiative cooling. Most of the Galactic 
volume (up to 70-80%) is filled with the hot (> 10^ K) low-density (~ 10"^ cm"^) 
ionized gas, interspersed by cold neutral and relatively dense clouds. Observationally it 
was difficult to determine reliable volume filling factors of the hot phase in our Galaxy 
due to the observational vantage point. However, it became evident that the small sur- 
face coverage of Hi holes in external galaxies (e.g.. Brinks & Bajaja|1986 1 argues for a 
much lower volume filling factor of the hot phase there. 



Furthermore, Population I stars typically bom in OB associations (see, e.g., Mc- 
Cray & Snow|1979| McCray & Kafatos|I987 l. Since more massive stars evolve more 
rapidly than low mass stars, the associations still exist when the first SNe occur. Hence, 
SNe are not randomly distributed over the entire galaxy. As a result of their motions, 
up to 40% of the stars explode in the field (see review by Ferriere|[2001 |l. Therefore, 



the amount of hot gas will mostly be concentrated in isolated pools (forming superbub- 



bles), which may burst through the thick disk (composed of Hi ( Lockman 1984 1 and Hii 
(Reynolds 1985) layers) into the halo like chimney funnels ( Norman & Ikeuchi|1989| ). 
The superbubble evolution in a static density stratified environment was modelled hy- 



drodynamically in two dimensions by, e.g., Tomisaka & Ikeuchi ( 1986 1, 
(^989), and |Tenorio-Tagle et al.| ( |T990l ). |Tbmi saka' ( 1 998|) and |S til eta 



Mac Low et al. 



(2009) devel- 



oped three-dimensional simulations of a superbubble evolving in a static magnetized 
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medium having a stratified structure along the direction perpendicular to the Galactic 
midplane. 

With the presence of such an extended thick disk, arguments were put forward 
that the break-out of bubbles and superbubbles could be inhibited (unless they occurred 
at a significant height above the disk) and only the most energetic superbubbles (SBs) 
would achieve blow-out of the disk (see Koo & McKee 1992) . Things might even 
become worse as superbubble (SB) break-out may be inhibited by a large-scale disk 
parallel magnetic field (e.g. Mineshige et al. 1993 Tomisaka|1998[ l. On the other hand, 
it seems paradoxical that tenestrial plasmas are so hard to confine, whereas ISM plas- 
mas should be magnetically tied down to the disk. In fact in both cases small and 
large scale instabilities destroy the highly symmetric configuration, like e.g. the Parker 
(1966) instability in the astrophysical context. Observationally, owing to the high sen- 
sitivity and large throughput of the ROSAT XRT ( |Trumper|[T983] l) and its PSPC in- 
strument, a number of normal spiral galaxies with soft X-ray halos were detected, e.g. 
NGC 891 ( |Bregman & Pildis|[l^94l ) and NGC 4631 (see [Vogler & Pietsch|[T996D . In 
some cases, even local correlations between Ho-, radio continuum and soft X-rays were 
found (Dettmar'1992 1, arguing for local outflows as it had been suggested in the Galac- 
tic fountain (Bregman|I980 tlKahn"1981), chimney CNorman & Ikeuchi"1989) and the 
Galactic wind model ( Breitschwerdt et al. 1991; Breitschwerdt & Schmutzler 1994). 

These models capture some of the structure but not all of the essential physics. 
Taking into account that the ISM is a turbulent and compressible system (von Weizsacker 
195 1), in which cooling and heating determine the physical state of the gas, more com- 



plex and sophisticated galactic disk (two-dimensional: jChiang & Prendergast 


1985[ 


Chiang & Bregman 1988 Rosen et al. 1993j) and disk-halo interaction two- ( 


Rosen 


& Bregman 1995) and three-dimensional (Korpi et al. 


1999a 


Avillez| 


2000[ J 


^villez 


& Breitschwerdt 2004 2005[ Joung & Mac Low 2006 


Melio 


i et al. 


2009 1 models 



were developed. In turn, some of these models have been used to follow the evolution 
of superbubbles within a realistic supernova-driven turbulent magnetized medium (in 
contrast to previous simulations that used a static medium) by [Korpi et al. ( 1999b I and 
Avillez & Breitschwerdt ( 2005| l (a field composed of a random and mean components 



with a total strength of 4.5 juG was used) finding that blow-out is much more likely than 
previously thought, mainly because bubbles evolve in an inhomogeneous background 
medium. However, Korpi et al.'s simulations are limited by the usage of a small grid 
extending only up to 1 kpc in the direction perpendicular to the disk in both sides of 
the midplane. Hence, the disk-halo-disk cycle could not be established nor followed. 
The simulations of Melioli et al. ( 2009 1 calculate the time evolution of Hi, Hii, Cii-Civ, 
and Oi-Oiu ions at T< 10" K and determine their contributions to the cooling function. 
For T> 10^ K they use a CIE cooling curve. Such setup has severe implications in the 
history and cooling of the plasma as will be discussed below. 

The evolution of superbubbles can be much improved by, instead of injecting en- 
ergy in a continuous way into a single point, identifying the missing stars in an associa- 
tion and follow in space and time all the stars during their main sequence life time until 
they explode. This methodology has been used in the disk-halo simulations, includ- 
ing those related to the evolution of the Local and Loop I superbubbles in a turbulent 
medium ( [Breitschwerdt & Avillez|[2006 ; Avillez & Breitschwerdt[[2009 ). In this way 
one can perform a detailed simulation of the properties of the bubbles and their spectro- 
scopic observables, allowing a direct comparison with observations. These simulations 
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have been successful in reproducing not only the spatial structure of the bubbles, but 
also the observed column densities of Li-hke ions Civ, Nv and Ovi and their ratios. 



3. Plasma Emission Modelling 



A further improvement in modelling the ISM comprises the full-blown non-equilibrium 
ionization structure (resulting from the ten most abundant elements), where the ioniza- 
tion, thermal and dynamical history of the plasma are fully nonlinearly coupled and 
tracked simultaneously both in space and time at high resolution. In order to achieve 
this, we developed a plasma emission code (hereafter EPEC - Eborae Plasma Emission 
Code) which can be coupled to any MHD software through the proper interfacing calls. 
EPEC is written in Fortran 2003 in an object-oriented way making large use of proce- 
dure pointers. It is prepared to run both on multi-core CPUs (using OpenMP) as well 
as on NVIDIA CPUs (graphics cards processor units) by means of CUDA Fortran. 



3.1. Abundances 



EPEC includes the ten most abundant elements in nature (H, He, C, N, O, Ne, Mg, Si, 



S and Fe) and the default solar abundances are those recommended by [Asplund et al. 
([2009. AGSS2009): logA(X///)o- -1.07 (He), -3.57 (C), -4.17 (N), -3.31 (O), -4.07 
(Ne), -4.40 (Mg), -4.49 (Si), -4.88 (S), and -4.50 (Fe). Other abundances (for com- 



parison studies with previously published results) are also available, e.g., Allen ( 1973 1, 
Anders & Grevesse| (fT989 | AG1989) and I Grevesse"etaLl (|2007 . GAS2007). The latter 
are used by Gnat & Sternberg ( 2007 ), but with the Ne overabundance (log A(Ne/H)o = 
-3.71) ofl Drake, &Tbsta| ( |2005) . GAS2007 and AGSS2009 propose C, N, O and Ne 
abundances smaller than those recommended by AG 1989, but AGSS2009 increases 
slightly those values from GAS2007. This variation in the abundances results from 
the improvement on three-dimensional hydrodynamical solar model atmospheres, that 
include a relaxation assumption on the local thermodynamic equilibrium and improve- 
ments in the atomic and molecular data (see discussion in, e.g., AGSS2009). 



3.2. Atomic and Cooling Processes 

The adopted physical processes in this work are the electron impact ionization, excitation- 
autoionization, radiative and dielectronic recombination, charge-exchange recombina- 
tion and ionization reactions, continuum and line emissions. Electron impact ionization 



rates fits are taken from [Mattioh et al] ( [2007, MAT2007) for all H, He, C, N, O, Ne 



ions, Mgi-Mgin, Sii-Siviu, Si-Sv, and Fei-Fexi; The experimental data of Fogle et al. 



( 2008 1 is used to fit the rates for Cni, Niv, and Ov. Data for the remaining ions is taken 
from'Maz zotta et aT (1998] MAZ1998). Excitation-autoionization rates are taken from 
MAT2007 for Civ, Nv, Ovi, Neviii, Siui-Siiv, Sm, Sv and Fexi. For the remaining ions 



where EA is important we follow Amaud & Rothenflug (1985 , AR85). 

Radiative recombination rates are fitted following ^Verner & Ferland | (|1996|) and 



Gu (2003 1, for low charge ions. The fits parameters for bare through Na-like ions are 
taken from Badnelf (|2006|) and for other ions we follow [DereetaD ([2009]). Dielectronic 



recombination rates fits coefficients for H-like through Mg-like ions are taken from 
2006]), ]Zatsarinny et aL] ([2003] ]2004] [2006 . and references therein), [CoIgaiTet 



Badnell 



al. (2003 



2004]),|Altun etal. !(2004, 2006. 2007 ), Mitnik & Bad nell] (]2004]) an d Bautista 



& Badnell] ( |2007] ) with updates for Svi ( |Orban et al.||2009| ), Nevii ( ]Orban et al.||2008f 7 
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Feviii-Fe ix ( [Schmidt et"aLl|20 08), Fex-Fexi ([Lestinsky et al.||2009| ), Fexiv ( [Schmidt et" 
al.|2006 l, Fexv ( Lukfc et al.|2005| , Fexxiii ( Savin et al.|2006 1. For the remaining ions 
MAZ1998 data is used. 

Charge-exchange recombination (CER) with Hi rates for Hen-Hem, Nn-Nv, Om, 
Ov, Neni-Nev, Mgra-Mgv, Siv, Sni-Sv and Feui-Fev are taken from'Kingdon & Ferland[ 
( [T9961), Cn-Cvn ([Suno & Kato|2006j, On ([Spirko et al ,2003j , Oiv (Wang et al. 2003| , 



Sim ( [Clarke et al.[[1998| ), and Siiv ( [Bruhns et al.|[2008[ ). Fits to the rates of CER with 
Hei for Nm-Nv, Om, Ov, Neiv-Nev, Civ-Cv, Mgiv-Mgv, Siiv-Siv, Siv and Feiv-Fev are 
taken from Astrophysics Charge-Transfer Databa se (jWang et al.|2002a and references 



therein). Oii ([Zhao et al.l[2005ll. Qiv (TWu et al.l[2009D , Nem dZhao et al.l|2006l ), Sm 



( |Zhao et al.|2005| , Sv ( [Wang et al.|2002b| ), and Fevi-Fexiv ( [Cadez et al.|2003[ ). Charge 



exchange ionization (CEI) with Hn rates for Mgu, Sii-Siii are taken from AR1985, Ci, 
Ni, Mgi, Si and Fei-Feii are taken from Kingdon & Ferland ( 1996[) a nd Oi from Spirko 



etal.[(|2003[ ). CEI with Hen comprised data for Oi ( [Zhao et al.|2004[ ) and Sin ( [Wang et 
al.|2002b[ ), Cn, Nn, Sim and Sii-Siii (AR1985). 



Cooling rates include free-free emission with the averaged Gaunt factor by Karzas 
& Latter (1961), radiative and dielectronic recombination ([Cox & Tucker 1969) , line 



° — I 1 1 — i r 

emission in the range 1 A -610 (|Penston|l970 Jura & Dalgamo 1972; Dalgarno, & 



McCray|1972t[Kato|1976t[Stem et al.|1978||Gaetz & Salpeter 1983 ; Mewe et al. 1985^ 1 



and two-photon emission. Spectra calculations include line and continuum emissions. 
The latter comprises free-free, free -bound to the ground and excited states, and two- 
photon - using ls-2s transitions in H and He-hke ions ( Tucker & Gould|1966 Gronen- 
[schild & Mewe|1978t[Mewe et al.|1986l . 



3.3. Equations 

The time evolution of the ions fractions, where ionization and recombinations of ions 
of nuclear charge Z occur between neighbouring ionization stages z - 1, z and z -i- 1, is 
given by 



dnz^ 



= Jz,z- l«Z,z-l«e - (i"z,z + f<Z,z)nZ,zne + %z+l«Z,z+l«e, 



(1) 



where "Rz,? ^z,z the rates of recombination and ionization from state (Z, z) to 
(Z, z - 1) and (Z, z + 1), respectively; ?iz,z and are the ion density of element Z with 
effective charge z and electron density, respectively. This constitutes a tridiagonal ma- 
trix if the charge exchange reactions are not included. However, if charge exchange 
reactions are included (as done in the present paper) new off-diagonals terms are intro- 
duced, becoming 



<^Z,z + "^Z,z + 



'Z,v 



(2) 



where the uppers indices stand for radiative, dielectronic and charge exchange pro- 
cesses, respectively, and 



T r^eii , /-tea 

- <-Z,z + <-Z,z 



"Z,z 



(3) 



with eii and ea standing for electron impact and excitation- autoionization processes, re- 
spectively, and - is the particle density of other ions involved in the charge exchange 
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reactions. In the EPEC version referred to here, we do neither consider photoioniza- 



tion nor ionization due to suprathermal electrons (see, e.g., Schmutzler & Tschamuter 
[T993^ - this is the subject of a forthcoming paper. For the 10 elements considered in the 
EPEC, the number of ordinary dilferential equations including the neutrals amounts to 
112. Conservation of species (atoms and ions) implies that 



nz = 

z=0 



£«z,z. (4) 



In addition the system of equations must account for the mass 



= ^ «z + ne (5) 



fltot 

Z 



and charge conservations 



z 



^ ^ ^ «ZzZ (6) 



Z z=l 

The system of equations is closed by the energy balance equation 

dU Pdn 

— — ^ -tieriHA (7) 

at n at 

where A is the cooling function, i.e., it represents the radiative losses per unit emission 
measure resulting from bremsstrahlung, radiative and dielectronic recombination, col- 
lisional ionization, line emission due to excitation, two-photon emission, and charge- 
exchange reactions. In these equations «// is the hydrogen density and U is the internal 
energy density of the system. The internal energy comprises the contributions from 
thermal motions of the particles, represented by Uth, and the potential energy associ- 
ated with the ionization stages, that is the energy stored in or delivered from the high 
ionization stages of chemical elements ( Schmutzler & Tschamuter||1993 1. Hence, the 



U ^Uth + Y^Tj ''z,z 2 

Z z=l 



internal energy density is given by 

z ( z-i ^ 

(8) 

z'=0 

where Iz,z' is the ionization potential of an ion with nuclear charge Z and effective 
charge z'. The thermal part of the internal energy is linked to the pressure of the system 
through the equation of state P = (y - l)Uth- 

As initial conditions for static plasma calculations we assume a fully ionized gas at 
10^ K. As the temperature decreases, recombination and ionization rates are calculated, 
following the simultaneous implicit calculation of the ionization fractions (including 
neutrals) and charge equation - implying the inversion of a 113x113 matrix. Next the 
radiative losses, emission spectra and internal energy are calculated. 



4. CoUisional ionization Equilibrium vs. Non-Equilibrium ionization 



Figure [T] shows the normalized cooling function of a static plasma (that is, with no 
dynamics included) that cooled from 10^ to 10^ K under CIE and NEI (isochorically) 
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conditions calculated with EPEC and using AGSS2009 abundances. Variations by fac- 
tors of 6 at lO'^-^ K are observed between the two efficiencies with the CIE cooUne 
dominating over the NEI case between lO^'-^ K and 10^ K (bottom panel of Figure [Ijl. 
These variations are a consequence of the cooling efficiency due to the different emis- 
sion processes (bremsstrahlung, radiative and dielectronic recombination, coUisional 
ionization, line excitation and two-photon emission) being larger under CIE than NEI 
conditions for all elements taken into account. Figure [2] compares the contributions of 
these processes to the cooling per element (e.g, C, Ne and Fe) under CIE (top row) and 
NEI (bottom row) conditions. 

For T< 10"* ' K the cooling efficiency at the same temperature under NEI (iso- 
choric) conditions is larger than in the CIE case as a result of the delayed recombination 
of the plasma. As recombination lags behind, single and double ionized species exist 
at lower temperatures (right column in Figure |3]l, something that does not occur un- 
der CIE, because ionization and recombination is synchronized, and ffierefore neutrals 
form at temperatures near 10^ K (left column in Figure[3|. 

While in CIE the ionization fractions (left panel Figure [3]l depend only on the 
temperature and are sharply peaked, in NEI these same fractions (right panel Figure [3]l 
depend on the dynamical and thermal history of the plasma. The higher ionization 
stages recombine to lower ones and eventually (when T~ 10^-^ K) only the lowest 
stages are abundant. However, the qualitative behaviour of all ionic stages is not the 
same. The highest ionic stages decrease continuously, that is, they always recombine 
to the next lowest stages, while the lowest stages increase continuously, that is the 
next highest stage recombines to them; the intermediate stages have two peaks result- 
ing from the recombination of the next highest stage, but dielectronic recombination 



■| Q' I — I I I II lll| 1 I I II lll| 1 1 I 11 lll| 1 11 II lii| 1 I I II ill| 1 I I II lll| 1 I I II lu 
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Figure 1. Top panel: Normalized CIE (daslied line) and NEI (solid line) cool- 
ing functions as function of temperature and calculated with AGSS2009 solar abun- 
dances. Bottom panel: CIE and NEI cooling functions ratio. 




Figure 2. Contributions to the CIE (top row) and NEI (bottom row) cooling due to 
bremsstrahlung (black dashed line), radiative and dielectronic recombination (black 
and green solid lines, respectively), collisional ionization (blue line), line excitation 
(red Une) and two-photon emission (brown line). 



rapidly depletes it, leading to the formation of the next lower ionization stage. When 
dielectronic recombination is no longer effective, recombination from the next highest 
stage increases the ion amount. As soon as the highest stage is depleted, the next stage 
recombines to the next lowest stage. 

These differences between CIE and NEI become quite noticeable in the emission 
spectra of the plasma at different temperatures (Figure |4]l. The spectra include line 
emission (cyan lines) and continuum (free-free (dashed black line), free -bound (green 
line) and two-photon (red line) emissions). With the decrease in temperature from 10^ 
to 10"*^ the CIE and NEI spectra become quite different as result of the free-bound 
emission dominating the spectra up to 500 A at low temperatures. At high tempera- 
tures, above 10^ K the differences between the ionization structure under CIE and NEI 
conditions are small, and therefore the spectra in these two cases are similar. They 
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Figure 3. Temperature variation of the ionization structure of C (top) and Ne (Bot- 
tom) under CIE (left column) and NEI (right column) conditions for a gas cooling 
from 10'' K. 



differ appreciably when the recombination lags behind, typically around 10^ K, with 
the continuum at low wavelengths (i.e., A < lOOA) becoming free-bound dominated. 
Note that in the CIE case there is very little («; 10~^° erg cm"^ s~^ A"^) emission with 
increasing range of wavelengths and decrease in temperature, e.g, at 10^ K: A < 73A, 
lO'^'^ K:A< 179A and lO'*^ K: A < 435A. Hence, no X -ray emission is expected in the 
latter two cases under CIE unlike in NEI, where due to delayed recombination, there is 
emission to be expected in these ranges of wavelengths and temperatures (right panel 
Figure |4]l. 

5. Signature of the Initial Conditions in the Turbulent ISM - NEI Modelling 

In a turbulent supernova-driven ISM the state of the plasma, and therefore its ionization 
structure, is determined by the heating and cooling as well as by the flow dynamics. 
Hence, it is expected that the ionization structure of the plasma varies from place to 
place leading to a multitude of cooling functions in the computational box. 

We carried out three-dimensional hydrodynamical disk-halo interaction simula- 
tions (in a patch of the Galaxy located at the solar radius with an area of 1 kpc^ parallel 
to the Galactic midplane, and extending to ±10 kpc perpendicular to it, similar to those 



described in Avillez & Breitschwerdt (2007 1 using a resolution of 0.5 pc (the highest 
so far used for large scale ISM evolution) corresponding to an effective grid with 2000^ 
cells per kpc^ boxes) and including (i) local self-gravity, (ii) heat conduction and (iii) 
time-dependent evolution of the ionization structure (of H, He, C, N, O, Ne, Mg, Si, 
S, and Fe) at each cell of the grid using EPEC described above; (iv) the revised solar 
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Figure 4. Spectra of a plasma that cooled from 10^ K under CIE (left column) and 
NEl (right column) conditions at 10^^ K (top panel) and 10^ K (bottom panel). The 
total continuum (black solid line) results from free-free (black dashed line), free- 
bound (green line) and two-photon (red line) emissions. Line emission is shown in 
cyan for the different wavelengths. Note the striking differences between CIE and 
NEl emission spectra. 



abundances by Asplund et al. (20091 are used; (v) Periodic boundary conditions are 
used along the vertical faces, while outflow boundary conditions are set at the top and 
bottom of the grid. 

The left panel of Figure[5]displays the density distribution in the Galactic midplane 
at evolution time 400 Myr, while the right panel shows the regions with temperatures 
between 10'^'^ K and 10^'' K. The displayed time is long enough for (i) the signature 
of the initial plasma conditions at f = to be wiped out (it takes some 80 Myr for 
this to happen; , Avillez & Breitschwerdt|2004] ), (ii) the disk-halo-disk circulation to be 
fully established, and (iii) the system to reach dynamical equilibrium in a statistical 
sense (Avillez 2000| Avillez & Breitschwerdt 2004). The figure shows in great detail 
the topology of the frothy and turbulent ISM permeated with hot gas and having low 
temperature filamentary structures. 

Bubbles and superbubbles dominate the landscape, but the volume filling factor 
(that is, the fraction of disk volume) of the hot gas is merely - 20% for the Galactic 
supernova rate. Dark blue regions have temperatures > 10^ K resulting from recent 
supernova activity, whereas darkest red regions are molecular clouds with an excess 
density of > 1000 cm~^ and temperature T< 100 K. Cold gas (red regions) is formed as 
a result of shock compressed layers and cooling instabilities in the flow. The details of 
these simulations due to the high resolution allow us to identify turbulent motions and 
their effects at small scales. Superbubbles and bubbles have turbulent flows of material 
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Figure 5. Left panel: Non-equilibrium ionization density distribution in the 
Galactic midplane at evolution time 400 Myr The color scale refers to the loga- 
rithm of the number density. Cold (high density) gas is represented by red while hot 
(low density) gas is shown in blue. Right panel: Midplane gas temperature between 
10^'^ and 10*' ' K at 400 Myr of evolution. The labels refer to a selection of sites 
having temperatures of 10^ K (labels A through E), 10^^ K (labels F through J) and 
10^ K (labels K through O). 



as can be seen in the cavities, which are crossed by sheet like structures with variable 
geometries. These turbulent flows are responsible for the redistribution of energy inside 
and outside of the cavities. 

Because the plasma keeps a record of its history, and therefore, of its initial condi- 
tion (temperature, total pressure and density - including mass fluxes from neighbouring 
locations), the ionization structure of the plasma varies from place to place leading to a 
multitude of cooling functions for the same initial temperature as can be seen in the left 
panel of Figure [6j which displays the cooling paths for gas having an initial tempera- 
ture of 10^ K and located at sites K through O, respectively (shown in the right panel of 
Figure [5]). The cooling efficiency fails to match that calculated under NET (isochoric) 
conditions (black dashed line) all the way down to 10^ K (left panel of Figure [6]l, be- 
cause the latter has less "potential energy" stored in high ionization stages. 

An important consequence of this variability seen in the ionization structure evo- 
lution and cooling paths is the occurrence of X-ray emission, through free-bound tran- 
sitions, at low temperatures, becoming larger than the corresponding NET emission in 
a static plasma. The NET spectra were calculated for a gas cooling in a time-dependent 
fashion from 10^ K. The right panel of Figure |6] shows that the free-bound emission 
from the K through O sites (red, green, blue, brown and magenta lines in the left panel) 
dominates the NET (solid black line) and CIE (dashed black line) static plasma emis- 
sion with decreasing temperature. This is a clear indication that recombination from 
the continuum is not following the cooling of the gas. 
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6. Discussion and Final Remarks 

In this paper we emphasize the importance in ISM simulations to follow the ionization 
structure of a plasma in a time-dependent fashion, coupled self-consistently to the dy- 
namics. In fact, the circulation of gas between the disk and halo is a dynamic process, 
which involves a time scale, that can be much shorter than any of the microphysical 
time scales due to ionization and recombination. The gas escaping into the halo has an 
initial temperature well in excess of 10^ K, where the assumption of coUisional ioniza- 
tion equilibrium (CIE) is approximately valid. As the hot plasma expands away from 
the disk it will cool adiabatically thereby reducing its temperature and density. The 
recombination of highly ionized species lags behind and occurs mainly at considerable 
heights from the disk. 

New state of the art simulations of the ISM, focusing on the detailed description 
of the plasma ionization structure in a supernova driven ISM show that: (i) in a dy- 
namic ISM, the ionization structure and, therefore, the cooling function, varies with 
time and from place to place, depending on the initial conditions and its history (a re- 
sult in accordance to previous discussions by |Kafatos| ( | 1 97 3| ), [Shapiro, & Moore| ( jl976| l 



and Sutherland & Dopita ( 1993)) regarding static plasmas, i.e., with no dynamics in- 



cluded), (ii) the cooling path can be quite different for gas even with the same initial 
temperature, but having different densities and pressures, (iii) this path may not follow 
the one predicted by the pure plasma emission calculations, that is, without the dynam- 
ics included, and (iv) as a consequence. X-ray emission occurs at temperatures < 10^ 
K. This is a consequence of the important emission contribution of delayed recombi- 




Figure 6. Left panel: Normalized cooling functions at sites K through O shown 
in the right panel of Figure [s] having an initial temperature of 10^ K. Right panels: 
Free-bound emission (color curves) from gas having an initial temperature of 10^ K, 
located at sites K-O and having the cooling curves displayed in the left panel. Solid 
and dashed black lines represent the emission expected from NEl and CIE static (i.e., 
no dynamics) plasma at the temperatures shown in each panel. 
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nation, arising from a severe mismatch of recombination and dynamical time scales of 
the plasma. Hence, as T decreases, the emissivity becomes much larger, by more than 
an order of magnitude, than that predicted by CIE at the same temperature. 

A fundamental consequence of the previous discussion, is that the cooling at any 
point in the ISM depends clearly on the history and dynamics of the plasma, that is, 
the ionization structure keeps a record of its initial conditions and the thermodynamic 
path it has taken. As the ionization structure varies from place to place there exists 
a multitude of different cooling functions in the ISM. Furthermore, with delayed re- 
combination having an important role in the cooling of overionized gas, i.e. gas which 
once had a high temperature during its history. X-ray emission at low temperatures is 



expected. Consequently, the mismatch of observed plasmas, e.g. by DXS (Sanders et 



[al.||2001; , XQC ( .McCammon et al.|[2002] ), EUVE ( |JeIinsky et a"L][T995^ , with standard 



CIE emission models (even for multi-temperature fits) may eventually result from not 
correctly describing the NEI structure of the plasmas. 

In addition, the amount of singly ionized ions which are important for molecular 
cloud chemistry, will also be affected by the timescales for recombination and cooling 
paths. This implies that a careful study of the ionization structure at higher temperatures 
is needed when one is dealing with molecular chemistry in the ISM. 
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